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I.  INTRODUCTION 


The  aerodynamics  of  a  family  of  flared  projectiles  has  been  investigated  experimentally^'  2 
in  the  Ballistic  Research  Laboratory^  Aerodynamics  Range.  These  projectiles  are  possible 
concepts  for  a  25mm  training  round.  Predictive  techniques  have  also  been  applied  to  examine 
the  drag  and  static  pitch-plane  aerodynamics  of  these  rounds. 2  In  this  report,  an  approach 
for  predicting  the  pitch-damping  aerodynamics  coefficients  has  been  applied  to  this  family 
of  projectiles.  This  approach  has  been  previously  applied  to  predict  the  pitch-damping  of 
both  finned  projectiles^  and  axisymmetric  shell.'* 

A  schematic  of  the  baseline  cone-cylinder-flare  projectile  configuration  is  shown  in  Fig¬ 
ure  1.  Each  of  the  projectiles  examined  here  has  the  same  cone-cylinder  forebody.  The 
forebody  has  a  slightly  truncated  conical  nose.  In  the  computations,  the  nose  is  modeled  as 
a  sharp  tipped  cone.  The  cylindrical  portion  of  the  body  also  has  a  number  of  subcaliber 
grooves  which  permit  the  launch  loads  to  be  transferred  from  the  sabot  to  the  projectile 
during  launch.  These  grooves  are  not  modeled  in  the  computations  presented  here. 

Various  afterbodies  have  been  analyzed  both  experimentally  and  computationally. 
Schematics  of  the  afterbodies  are  shown  in  Figure  2.  Additional  dimensions  are  listed  in 
Table  1.  The  configurations  CS-V4-2  through  CS-V4-5  have  a  1-caliber  afterbody  extention 
added  to  the  b<iseline  configuration,  CS-V4-1.  The  angle  of  inclination  of  the  conical  exten- 
tions  for  configurations  CS-V4-2  through  CS-V4-5  are,  respectively,  6°  (simple  extention  of 
original  flare),  0°  (cylindrical  skirt),  12®  (steeper  flare),  and  —6°  (boattail).  Configuration 
CS-V4-6  consists  of  9.37®  flare  which  has  been  machined  to  produce  a  square  cross  section 
over  the  last  caliber  of  the  body.  Configuration  CS-V4-7  is  identical  to  the  boattailed  con¬ 
figuration  CS-V4-5,  except  that  four  12®  fins  have  been  added  to  the  boattailed  portion  of 
the  body.  The  fins  are  0.153  calibers  thick.  The  final  configuration,  CS-V4-8,  is  identical  to 
the  baseline  configuration,  except  that  four  boundetry  layer  strakes  have  been  added  to  the 
flared  portion  of  the  body.  The  strakes  are  0.153  calibers  in  height  and  width. 

The  flared  portion  of  the  bodies  have  been  hollowed  to  shift  the  center  of  gravity  forward 
to  produce  statically  stable  configurations.  There  is  some  evidence  that  the  hollow  base  may 
have  some  effect  on  the  aerodynamics.  In  this  investigation,  the  flow  field  in  the  base  region 
has  not  been  computed,  hence,  the  effect  of  the  hollow  base  on  the  aerodynamics  cannot  be 
assessed. 


II.  THEORETICAL  BACKGROUND 


By  applying  linear  flight  mechanics  theory  such  as  that  developed  by  Murphy,®  it  can 
be  shown  that  aerodynamic  side  force  and  moment  coefficients  acting  on  a  projectile  in 
steady  coning  motion  can  be  related  to  the  pitch-damping  force  and  moment  coefficients. 
Steady  coning  motion  is  defined  as  the  motion  performed  by  a  missile  flying  at  a  constant 
angle  with  respect  to  the  free-stream  velocity  vector  and  undergoing  a  rotation  at  a  constant 

•The  U.S.  Army  Ballistic  R^arch  Laboratory  was  deactivated  on  30  September  1992  and  subsequently 
became  a  part  of  the  U.S.  Army  Research  Laboratory  (ARL)  on  1  October  1992. 
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angular  velocity  about  a  line  parallel  to  the  free-stre<im  velocity  vector  and  coincident  with 
the  projectile  center  of  gravity  (Figure  3).  Coning  motion  is,  in  fact,  a  specific  combination 
of  two  orthogonal  planar  pitching  motions.  The  use  of  steady  coning  motion  to  determine  the 
pitch-damping  aerodynamic  coefficients  provides  an  interesting  and  cost-effective  approach 
for  determining  the  aerodynamics  which  are  normally  tissociated  with  unsteady  or  time- 
dependent  motions. 

Previously,  Tobak,  Schiff,  and  Peterson®  examined  the  aerodynamics  of  bodies  of  revo¬ 
lution  in  coning  motion  and  proposed  that  the  nonlinear  aerodynamic  forces  and  moments 
acting  on  a  body  performing  large  amplitude  nonplanar  motions  could  be  composed  of  four 
characteristic  motions:  (1)  steady  angle  of  attack,  (2)  pitching  motion,  (3)  roiling  motion, 
and  (4)  coning  motion.  Typically,  the  linear  aerodynamic  force  and  moment  formulation  con¬ 
siders  only  forces  and  moments  due  to  the  first  three  motions,  and  assumes  that  a  nonplanar 
motion  can  be  described  by  the  vector  sum  of  two  independent  planar  motions.  The  addition 
of  coning  motion  allows  for  coupling  between  planar  motions  in  the  nonlinear  formulation. 
Their  nonlinear  theory  also  confirms  the  linear  theory  result  that  the  side  force  and  moment 
due  to  coning  motion  is  related  to  the  linear  pitch-damping  coefficients.  Schiff  and  Tobak^ 
performed  wind  tunnel  experiments  on  a  conical  body  undergoing  separate  or  combined 
spinning  and  coning  motions  to  provide  additional  validation  for  the  theory. 

As  part  of  the  development,  Schiff®  computed  the  supersonic  inviscid  flow  about  a  coni¬ 
cal  body  undergoing  coning  motion.  To  compute  the  flow  around  the  body  in  coning  motion, 
Schiff  made  use  of  a  rotating  coordinate  fr<ime.  Within  the  rotating  coordinate  frame  the 
flow  was  steady,  thus  the  steady  Euler  equations  could  be  solved.  The  governing  equations 
were  modified  to  include  the  centrifugal  and  Coriolis  force  terms.  His  computed  results  com¬ 
pared  well  with  experimental  results  and  with  estimates  of  pitch-damping  coefficients  using 
a  linear  theory.  Later  studies  by  Agarwal  and  Rakich®  emd  Lin^^also  employed  rotating 
coordinate  frames  to  compute  the  supersonic  viscous  flow  about  conical  bodies  in  coning 
motion.  More  recently,  Weinacht  and  Sturek^  performed  computations  for  finned  projectiles 
in  coning  motion  to  determine  the  pitch-damping  coefficients. 

In  each  of  these  efforts,  the  pitch-damping  coefficients  were  determined  from  the  side 
moment  due  to  lunar  coning  motion.  Lunar  coning  motion  is  a  specific  form  of  coning  motion 
in  which  the  body  does  not  rotate  with  respect  to  the  pitch  plane.  The  relation  between 
the  side  moment  due  to  lunar  coning  motion  and  the  pitch-damping  moment  is  shown  in 
equation  1  without  derivation.  This  expression  was  obtained  using  linear  flight  mechanics 
theory  (details  can  be  found  in  reference  [4]  )  and  is  essentially  the  same  as  that  obtained 
earlier  by  Schiff  and  Tobak'^  by  applying  their  nonlinear  theory  to  bodies  of  revolution. 

Cn  =  <(||)(7Cn„  +  (C„,+7C„J)  (1) 

The  side  moment  in  the  coning  frame  of  reference,  Cn,  varies  linearly  with  the  coning  rate, 
<f>,  and  the  sine  of  the  angle  of  attack,  S,  and  is  a  function  of  the  pitch-damping  moment 
coefficient,  Cm,  +  'fCma,  and  the  Magnus  moment  coefficient,  Cn^a-  The  cosine  of  the  angle 
of  attack,  7,  is  often  dropped  because  it  is  close  to  one  in  the  region  where  linear  theory  is 
valid.  In  many  cases,  the  Magnus  moment  coefficient  is  much  smaller  than  the  pitch-damping 
coefficient,  and  the  Magnus  moment  can  be  dropped  allowing  the  side  moment  due  to  lunar 
coning  to  be  directly  related  to  the  pitch-damping  moment  coefficient. 
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Recently,  Wei^  xht,  Sturek,  and  SchifF"*  computed  the  pitch-damping  for  axisymmetric 
shell  using  a  specific  combination  of  spinning  and  coning  motion  which  allows  the  pitch¬ 
damping  force  and  moment  coefficients  to  be  directly  related  to  the  side  force  and  moment 
due  to  this  motion,  as  shown  in  equation  2. 

Cn  =  <(|J)|C„.  +  -,C„.l  (2) 

In  the  coning  coordinate  frame,  the  body  will  spin  in  the  opposite  direction  of  the  coning 
motion,  where  the  spin  rate  in  the  coning  frame,  pc}  is  related  to  the  coning  rate  in  the 
following  manner;  pcj  =  — Despite  the  simplicity  of  this  expression,  the  Magnus  problem 
has  not  been  entirely  removed  from  the  problem,  however.  This  is  because  the  motion  in 
the  coning  coordinate  frame  involves  coning  and  spinning  motions.  Thus,  any  approach, 
whether  it  be  computational  or  experimental,  which  uses  this  motion  must  be  capable  of 
modeling  both  of  these  effects. 

Similar  expressions  which  relate  the  side  force  due  to  coning  to  the  pitch-damping  force 
and  Magnus  force  can  be  derived  for  the  cases  of  lunar  coning  and  combined  spinning  and 
coning  motion. 


III.  COMPUTATIONAL  APPROACH 


Computation  of  the  viscous  flow  field  about  the  flared  projectile  configurations  was  ac¬ 
complished  by  solving  the  thin-layer  Navier-Stokes  equations  using  the  parabolized  Navier- 
Stokes  technique  of  Schiff  zuid  Steger.^* Using  the  parabolized  Navier-Stokes  technique,  com¬ 
putational  results  are  obtained  by  marching  through  the  grid  from  the  projectile  nose  to  the 
base.  This  technique  is  applicable  in  the  supersonic  flow  regime  and  requires  that  the  flow 
field  contain  no  regions  of  flow  separation  in  the  axial  direction.  Because  the  computational 
approach  requires  only  a  single  sweep  through  the  computational  grid,  it  is  very  efficient 
compared  with  time-marching  approaches  which  require  many  sweeps  through  the  grid. 


The  flow  field  predictions  of  the  projectile  in  steady  coning  motion  have  been  performed 
using  a  rotating  coordinate  frame  which  rotates  at  the  coning  rate  of  the  projectile.  The 
fluid  flew  relative  to  the  rotating  coordinate  frame  does  not  vary  with  time,  allowing  the 
steady  (non-time  varying)  Navier-Stokes  equations  to  be  applied.  To  implement  the  rotat¬ 
ing  coordinate  frame,  the  governing  equations  have  been  modified  to  include  the  effect  of 
centrifugal  and  Coriolis  forces.  The  steady  thin-layer  Navier-Stokes  equations  are  shown  in 
equation  3. 

dE  dF  dG  -  1  dS 

di^  ~  Re  dC, 


tk  A  A 

Here,  E,  F,  and  G  are  the  inviscid  flux  vectors,  5  is  the  viscous  flux  vector,  and  H  is  the 
source  term  containing  the  Coriolis  and  centrifugal  force  terms  which  result  from  the  rotating 
coordinate  frame.  Each  of  these  matrices  are  functions  of  the  dependent  variables  represented 
by  the  vector  q{p,  pu,  pv,  pw,  e)^,  where  p  and  e  are  the  density  and  the  total  energy  per  unit 
volume,  and  u,  v,  and  w,  are  the  velocity  components  in  x,  y,  and  z  directions.  The  pitch 
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plane  is  oriented  in  the  x  —  z  plane.  The  flux  terms  are  shown  in  equation  4. 


^=7 


where 
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law. 


The  pressure,  p,  can  be  related  to  the  dependent  variables  by  applying  the  ideal  gas 


p  =  (7  -  1)  e  -  ^9^ 


(11) 
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The  turbulent  viscosity,  /*<,  which  appears  in  the  viscous  matrices,  was  computed  using  the 
Baldwin- Lomax  turbulence  model. 

The  thin-layer  equations  are  solved  using  the  parabolized  Navier-Stokes  technique  of 
SchifF  and  Steger.^’ Following  the  approach  of  SchifF  and  Steger,  the  governing  equations, 
which  have  been  modified  here  to  include  the  Coriolis  and  centrifugal  force  terms,  are  solved 
using  a  conservative,  approximately  factored,  implicit  finite-difference  numerical  algorithm 
as  formulated  by  Beam  and  Warming. ‘^Further  details  on  the  implementation  of  the  source 
terms  due  to  rotating  coordinate  frame  can  be  found  in  references  [3]  and  [4]. 

The  computations  presented  here  were  performed  using  a  shock-fitting  procedure  re¬ 
ported  by  Rai  and  Chaussee.^'^This  procedure  solves  the  five  Rankine-Hugoniot  jump  con¬ 
ditions,  two  geometric  shock-propagation  conditions,  and  one  compatibility  equation  to  de¬ 
termine  the  values  of  the  five  dependent  variables  immediately  behind  the  shock,  ais  well 
as  the  position  of  the  shock.  By  including  the  implicit  part  of  the  source  term  due  to  the 
rotating  coordinate  frame  in  the  circumferential  inversion,  the  shock-fitting  procedure  of  Rai 
and  Chaussee  can  be  used  without  modification,  £is  long  as  the  correct  free-stream  conditions 
are  specified  as  shown  in  equation  12  in  nondimensional  form. 

9  =  1 

pu  =  A/ooCoso: -f- yfFcSina 
pv  =  ftc('2cosa  —  (z  —  a:cg)sin«) 
pw  =  Moosina  —  yVlc  cos  a 

e  =  poo/(7  -  1)  +  ^{(A^oocosa-I- yficsina)^ 

-f(  17c (2: COSO  —  (a:  —  Xcjjsina))^  -I-  {M^sina  —  j/H,  cosq)^}  (12) 

On  the  body  surface,  constant  (ambient)  wall  temperature,  no-slip  boundary  conditions  were 
applied.  For  the  cases  involving  spin,  the  horizontal  and  vertical  velocities  at  the  projectile 
surface  were  appropriately  modified  to  account  for  the  tangential  velocity  due  to  solid  body 
rotation. 


IV.  RESULTS 


Using  side  force  and  moment  due  to  steady  lunar  coning  motion,  predictions  of  the  pitch¬ 
damping  force  and  moment  coefficients  have  been  made  for  the  family  of  flare  projectiles. 
The  computations  were  performed  over  a  range  of  Mach  numbers  (A/00  =  3.0  to  4.5)  for 
free-flight  (sea-level)  atmospheric  conditions. 

The  computed  variation  of  the  side  moment  coefficient  with  coning  rate  for  configuration 
CS-V4-1  at  Mach  4  and  2®  angle  of  attack  is  shown  in  Figure  4.  The  variation  is  linear  over 
the  range  of  coning  rates  examined  here,  and  the  slope,  Cn^,  can  be  determined  without  any 
ambiguity.  The  range  of  coning  rates  considered  here  encompasses  the  pitching  frequency  of 
the  projectile. 

Figure  5  shows  the  variation  of  the  side  moment  slope,  Cn^,  with  aagle  of  attack  for 
configuration  CS-V4-1.  The  angle  of  attack  variation  is  seen  to  be  linear  up  to  about  5°  angle 
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of  attack.  At  angles  of  attack  less  than  5°,  the  normalized  side  moment  slope,  varies  by 

only  a  few  percent  as  the  angle  of  attack  changes.  Beyond  5°  angle  of  attack,  the  side  moment 
shows  a  strong  nonlinear  behavior  with  angle  of  attack.  In  the  development  of  the  relation 
between  the  side  moment  due  to  coning  and  the  pitch-damping  moment,  the  aerodynamic 
coefficients  were  assumed  to  vary  linearly  with  angle  of  attack.  At  angles  of  attack  greater 
than  5®,  the  previously  presented  relation  between  the  pitch-damping  coefficients  and  the 
side  force  £ind  moment  due  to  coning  motion  may  no  longer  be  valid. 

From  the  side  moment  slope  normalized  by  the  sine  of  the  angle  of  attack,  the  pitch¬ 
damping  coefficients  as  a  function  of  Mach  number  were  determined  for  each  of  the  eight 
flared  configurations.  The  predictions  were  obtained  at  5®  angle  of  attack.  Figures  6-13  dis¬ 
play  the  predictions  of  the  pitch-damping  moment  coefficients,  as  well  as  the  pitch-damping 
moment  coefficients  determined  from  the  free  flight  range  tests.  The  predictions  of  the 
pitch-damping  moment  are  within  the  accuracy  of  the  range  measurements.  The  predicted 
pitch-damping  moment  coefficients  for  each  of  the  configurations  show  a  weak  dependence 
on  Mach  number  over  the  range  of  Mach  numbers  considered  here. 

Figure  14  shows  a  comparison  of  the  pitch-damping  moment  coefficients  for  each  of  the 
eight  configurations  at  Mach  4.  Both  the  PNS  predictions  and  the  range  data  are  shown. 
Each  of  the  bodies  with  the  conical  extentions,  configurations  CS-V4-2  through  CS-V4-5, 
have  larger  pitch-damping  coefficients  compared  with  the  baseline  configuration  CS-V4-1. 
Computational  predictions  show  a  consistent  increase  in  the  pitch-damping  for  the  bodies 
with  the  conical  extentions  with  the  boattailed  configuration  having  lowest  pitch-damping 
coefficient  and  the  steepest  flare  having  the  highest  pitch-damping  coefficient.  These  trends 
are,  for  the  most  part,  reflected  by  the  range  data. 

The  finned  configuration,  which  is  identical  to  the  boattailed  configuration  except  that 
four  12®  swept  fins  have  been  added  to  the  boattailed  portion  of  the  body,  shows  a  modest 
increase  in  the  damping  over  the  boattailed  configuration.  The  finned  configuration,  however, 
produces  significantly  less  pitch-damping  than  the  configuration  with  the  12®  flare  extention. 
Again,  these  trends  are  reflected  by  the  r<inge  data.  The  configuration  with  the  boundary 
layer  strakes  also  produces  a  modest  increase  in  the  damping  compared  with  the  baseline 
configuration  which  has  no  strakes.  The  square  base  configuration  (CS-V4-6),  which  has 
the  same  base  area  as  the  beiseline  configuration  (CS-V4-1)  and  the  configuration  with  the 
cylindrical  skirt  (CS-V4-3),  produces  more  damping  than  the  baseline  configuration  and 
slightly  more  damping  than  the  cylindrical  skirt. 

In  addition  to  pitch-damping  moment  coefficient,  predictions  of  the  the  pitch-damping 
force  coefficient  have  been  made.  The  pitch-damping  force  coefficient  can  be  determined 
from  the  side  force  variation  with  coning  rate  normalized  by  the  sine  of  the  angle  of  attaek. 
The  predicted  variation  of  the  pitch-damping  force  coefficient  with  Mach  number  for  config¬ 
uration  CS-V4-1  is  shown  in  Figure  15.  This  coefficient  is  seen  to  decrease  with  increasing 
Mach  number  for  the  range  of  Mach  numbers  considered  here.  The  Mach  number  variation 
was  typical  of  that  observed  for  each  of  the  eight  configurations.  The  pitch-damping  force 
coefficient  produces  little  effect  on  the  motion  of  the  projectile  and  is  very  difficult  to  deter¬ 
mine  from  the  in-flight  projectile  motion.  For  this  reason  no  range  data  are  shown.  Although 
the  pitch-damping  force  coefficient  is  not  required  to  determine  the  in-flight  motion,  this  co¬ 
efficient  can  be  useful  in  determining  the  variation  of  the  pitch-damping  moment  coefficient 
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with  change  in  center  of  gravity.  For  this  reason,  a  predictive  approach  for  determining  this 
coefficient  is  of  interest. 

Determining  the  pitch-damping  coefficients  from  the  side  force  and  moment  due  to  lu¬ 
nar  coning  motion  requires  that  the  Magnus  force  and  moment  be  determined  from  another 
source  or  neglected.  For  the  ajcisymmetric  configurations  (CS-V4-1  to  CS-V4-5),  the  Magnus 
force  and  moment  have  been  predicted  using  the  PNS  approach.  These  calculations  were 
performed  with  the  body  spinning  at  angle  of  attack  and  in  the  absence  of  coning  motion. 
Magnus  predictions  for  the  non-axisymmetric  geometries  could  not  be  made  because  the 
combination  of  spin  and  angle  of  attack  for  these  bodies  produces  a  time-dependent  flow 
field.  Figures  16  and  17  shows  the  predicted  Magnus  force  and  moment  as  a  function  of 
Mach  number  for  each  of  the  axisymmetric  configurations.  The  configuration  with  the  boat- 
tailed  afterbody  shows  the  largest  positive  Magnus  moment,  while  the  12°  flared  afterbody 
shows  the  largest  negative  Magnus  moment.  These  two  configurations  also  produce  the 
smallest  and  largest  Magnus  force  coefficients.  Figure  18  shows  a  comparison  of  the  Magnus 
moment  coefficients  by  configuration  at  Mach  4.  Range  results  are  shown  for  each  of  the 
eight  configurations,  while  computational  results  are  shown  for  the  axisymmetric  configura¬ 
tions.  The  computational  data  are  bracketed  by  the  range  data  for  each  of  the  axisymmetric 
configurations,  and  similar  trends  are  shown  by  the  computational  and  experimental  results. 
Both  the  computational  and  experimental  results  reveal  that  the  Magnus  moment  is  small  in 
comparison  to  the  pitch-damping  moment  coefficient  for  the  configurations  examined  here. 
The  computational  predictions  also  confirm  that  the  Magnus  force  is  small  in  relation  to  the 
pitch-damping  force  coefficient.  This  result  demonstrates  that,  by  ignoring  the  contribution 
from  the  Magnus  coeflScients,  the  pitch-damping  coefficients  can  be  determined  directly  from 
the  side  force  and  moment  due  to  lunar  coning  motion  with  little  effect  on  the  accuracy  of 
the  prediction. 

For  the  axisymmetric  configurations,  calculations  were  also  performed  using  combined 
spinning  and  coning  motion.  The  predicted  pitch-damping  coefficients,  which  are  obtained 
directly  from  the  normalized  side  force  and  moment  (Equation  2),  showed  differences  of  less 
than  0.02  %  with  results  obtained  by  subtracting  the  Magnus  coefficients  from  the  normalized 
side  force  and  moment  coefficients  due  to  lunar  coning  motion  (Equation  1).  This  result 
demonstrates  the  lack  of  coupling  between  coning  motion  and  spin  over  the  range  of  spin 
and  coning  rates  considered  here  for  small  angles  of  attack. 

The  Magnus  moment  predictions  shown  here  are  interesting  because  the  predictions 
show  that  the  sign  of  the  Magnus  moment  is  dependent  on  the  projectile  geometry.  The 
Magnus  moment  results  from  an  interaction  between  spin  and  angle  of  attack  which  causes  a 
three  dimensional  boundary  layer  displacement  effect.  The  magnitude  of  the  effect  is  greatest 
near  the  aft  end  of  the  body  where  the  boundary  layer  is  thickest.  Using  the  ballistic  force- 
moment  sign  conventions,  the  result  for  axisymmetric  shell  is  typically  a  negative  Magnus 
force  and  a  positive  Magnus  moment  for  cylindrical  afterbodies  with  or  without  a  boattail  in 
supersonic  flow.  For  three  of  the  flared  afterbodies,  a  negative  Magnus  moment  is  predicted, 
with  the  reversal  in  Magnus  moment  occurring  over  the  flared  portion  of  the  body.  A  possible 
explanation  for  this  phenomenon  is  that  the  boundary  layer  zisymmetry  produced  by  angle 
of  attack  and  spin  causes  less  side  force  on  the  side  of  the  flare  where  the  boundary  layer  is 
thickest  due  to  a  blanketing  effect.  The  result  is  a  positive  contribution  to  the  Magnus  force 
and  because  this  occurs  at  the  aft  end  of  the  body,  a  negative  contribution  to  the  Magnus 
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moment  is  produced.  This  explanation  appears  to  be  supported  by  calculations  which  show 
that  for  these  flared  bodies,  lift  and  pitching  moment  increcises  as  the  Reynolds  number 
decreases,  presumably  due  to  a  reduction  in  pressure  on  the  lee  side  of  the  body  caused  by 
the  thicker  boundary  layer. 


V.  CONCLUSION 


The  pitch-damping  coefficients  for  a  family  of  flared  projectiles  have  been  successfully 
computed  using  a  parabolized  Navier-Stokes  approach.  The  pitch-damping  force  and  moment 
coefficients  have  been  predicted  from  the  side  force  and  moment  due  to  steady  lunar  coning 
motion.  The  side  moment  is  seen  to  vary  linearly  with  coning  rate  across  the  range  of 
interest.  The  side  moment  also  varies  linearly  with  angle  of  attack  up  to  about  5°.  These 
predictions  provide  a  constraint  on  the  range  of  angles  of  attack  where  the  pitch-damping 
coefficients  can  be  determined  from  the  side  force  and  moment  due  to  coning.  The  pitch¬ 
damping  moment  coefficient  predictions  for  the  family  of  flared  projectiles  has  been  compared 
with  data  obtained  from  free-flight  aerodynamics  range  tests  and  the  comparisons  show  good 
agreement.  The  predictions  show  that  the  pitch-damping  moment  coefficient  for  these  flared 
projectiles  is  relatively  insensitive  to  Mach  number  at  the  range  of  supersonic  velocities 
examined  here.  Predictions  of  the  pitch-damping  force  coefficient  have  also  been  made,  but 
no  comparisons  with  range  data  are  shown  due  to  the  small  effect  of  this  coefficient  on  the 
inflight  motion  of  the  projectile.  The  Magnus  coefficients  determined  from  the  predictive 
approach  and  from  range  tests  is  seen  to  be  small  in  relation  to  the  pitch-damping  coefficients. 
By  ignoring  the  Magnus  coefficients,  the  pitch-damping  coefficients  can  be  determined  from 
the  side  force  and  moment  due  to  lunar  coning  motion  with  little  effect  on  the  accuracy  of 
the  prediction. 
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Table  1.  Physical  Dimensions  of  the  Family  of  Flared  Projectiles 


Configuration 

Flare 

Angle 

Flare 

Fineness 

Ratio 

Extention 
(1  Caliber) 

CG  Position 
(Calibers 
from  Base) 

Total 

L/D 

CS-V4-1 

3.51 

None 

7.55 

15.36 

CS-V4-2 

4.49 

None 

8.24 

16.34 

CS-V4-3 

wmm 

4.49 

6.0®  Flare 

8.36 

16.34 

4.49 

Cylindrical  Skirt 

8.13 

16.34 

CS-V4-5 

4.49 

12.0®  Flare 

8.35 

16.34 

CS-V4-6 

4.49 

6.0®  Boattail 

8.05 

16.34 

CS-V4-7 

ITiftTiai 

4.50 

Boattail  with  Fins 

8.36 

16.35 

CS-V4-8 

3.50 

Strakes 

7.52 

15.35 
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ALL  DIMENSIONS  IN  CALIBERS  (ONE  CALIBER  =  8.28 
Figure  1.  Schematic  of  baseline  projectile  configuration 


CS-V4-1  CS-V4-2  CS-V4-3 


CS-V4-4  CS-V4-6 


Figure  2.  Schematic  of  projectile  afterbodies 


Figure  3.  Schematic  of  coning  motion 


MACH  NUMBER 


Figure  6.  Variation  of  the  pitch-damping  moment  coefficient  with  Mach  number,  CS-V4-1 


MACH  NUMBER 


Figure  7.  Variation  of  the  pitch-damping  moment  coefficient  with  Mach  number,  CS-V4-2 
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MACH  NUMBER 


Figure  8.  Variation  of  the  pitch-damping  moment  coefficient  with  Mach  number,  CS-V4-3 


MACH  NUMBER 


Figure  9.  Variation  of  the  pitch-damping  moment  coefl&cient  with  Mach  number,  CS-V4-4 
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MACH  NUMBER 

Figure  10.  Variation  of  the  pitch- damping  moment  coefficient  with  Mach  number,  CS-V4-5 


CS-V4-5 


Figure  11.  Variation  of  the  pitch-damping  moment  coefficient  with  Mach  number,  CS-V4-6 


MACH  NUMBER 


Figure  12.  Variation  of  the  pitch-damping  moment  coefficient  with  Mach  number,  CS-V4-7 


MACH  NUMBER 


Figure  13.  Variation  of  the  pitch-damping  moment  coefficient  with  Mach  number,  CS-V4-8 
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CONFIGURATION 


Figure  14.  Pitch-damping  moment  coefficient  versus  configuration,  Mach  4 


MACH  NUMBER 


Figure  15.  Variation  of  the  pitch-damping  force  coefficient  with  Mach  number,  CS-V4-1 
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MACH  NUMBER 


Figure  16.  Variation  of  the  Magnus  force  coefficient  with  Mach  number,  CS-V4-1  to 
CS-V4-5 


MACH  NUMBER 

Figure  17.  Variation  of  the  Magnus  moment  coefficient  with  Mach  number,  CS-V4-1  to 
CS-V4-5 


17 


0  range  data 


a  cn  cm  C3]  cTD  CT  ca  S 


REFERENCES 


1.  Celmins,  I.,  “Aerodynamic  Characteristics  of  Fin  and  Flare-Stabilized  25mm  XM910 
Prototypes,”  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD, 
Technical  Report  BRL-TR-2882,  December  1987.  (AD  191683) 

2.  Danberg,  J.E.,  Sigal,  A.,  and  Celmins,  I.,  “Aerodynamic  Characteristics  of  a  Family 
of  Cone-Cylinder-Flare  Projectiles,”  Journal  of  Spacecraft  and  Rockets,  Vol.  27,  No.  4, 
July-August  1990,  pp  335-360.  Also  see,  Danberg,  J.E.,  Sigal,  A.,  and  Celmins,  I.,  “Pre¬ 
diction  and  Comparison  with  Measurements  of  the  Aerodynamic  Characteristics  of  Flaxe- 
Stabilized  XM910  Prototypes,”  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Prov¬ 
ing  Ground,  MD,  Technical  Report  BRL-MR-3752,  May  1989.  (AD  208099) 

3.  Weinacht,  P.,  and  Sturek,  W.B.,  “Navier-Stokes  Predictions  of  Pitch  Damping  for  Finned 
Projectiles  Using  Steady  Coning  Motion,”  AIAA  Paper  90-3088,  AIAA  8th  Applied  Aero¬ 
dynamics  Conference,  Portland,  OR,  August  20-22,  1990. 

4.  Weinacht,  P.,  Sturek,  W.B.,  and  Schiff,  L.B.,  “Navier-Stokes  Predictions  of  Pitch  Damp¬ 
ing  for  Axisymmetric  Shell  Using  Steady  Coning  Motion,”  AIAA  Paper  91-2855,  AIAA 
Atmospheric  Flight  Mechanics  Conference,  New  Orleans,  LA,  August  12-14,  1991. 

5.  Murphy,  C.H.,  “Free  Flight  Motion  of  Symmetric  Missiles,”  U.S.  Army  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD,  Report  No.  1216,  July  1963.  (AD  A442757) 

6.  Tobak,  M.,  Schiff,  L.B.,  and  Peterson,  V.L.,  “Aerodynamics  of  Bodies  of  Revolution  in 
Coning  Motion,”  AIAA  Journal.  Vol.  7,  No.  1,  January  196G,  pp.  95-99. 

7.  Schiff,  L.B.,  and  Tobak,  M.,  “Results  from  a  New  Wind-Tunnel  Apparatus  for  Studying 
Coning  and  Spinning  Motions  of  Bodies  of  Revolution,”  AIAA  Journal.  Vol.  8,  No.  11, 
November  1970,  pp.  1953-1957. 

8.  Schiff,  L.B.,  “Nonlinear  Aerodynamics  of  Bodies  in  Coning  Motion,”  AIAA  Journal.  Vol. 
10,  No.  11,  November  1972,  pp.  1517-1522. 

9.  Agarwal,  R.,  and  Rakich,  J.V.,  “Computation  of  Supersonic  Laminar  Viscous  Flow  Past  a 
Pointed  Cone  at  Angle  of  Attack  in  Spinning  and  Coning  Motion,”  AIAA  Paper  78-1211, 
AIAA  11th  Fluid  and  Plasma  Dynamics  Conference,  Seattle,  WA,  July  1978. 

10.  Lin,  T.C.,  “A  Numerical  Study  of  the  Aerodynamics  of  a  Reentry  Vehicle  in  Steady 
Coning  Motion,”  AIAA  Paper  78-1358,  AIAA  Atmospheric  Flight  Mechanics  Conference, 
Palo  Alto,  CA,  August  1978. 

11.  Schiff,  L.B.,  and  Steger,  J.L.,  “Numerical  Simulation  of  Steady  Supersonic  Viscous 
Flow,”  AIAA  Journal.  Vol.  18,  No.  12,  December  1980,  pp.  1421-1430. 

12.  Baldwin,  B.S.,  and  Lomax,  H.,  “Thin  Layer  Approximation  and  Algebraic  Model  for 
Separated  Turbulent  Flows,”  AIAA  Pai)er  78-257,  16th  Aerospace  Sciences  Meeting, 
January  1978. 

13.  Beam,  R.,  and  Warming,  R.F.,  “An  Implicit  Factored  Scheme  for  the  Compressible 
Navier-Stokes  Liquations,”  AIAA  Journal.  Vol.  16,  No.  4,  1978,  pp.  85-129. 


14.  Rai,  M.M.,  and  Chaussee,  D.S.,  “New  Implicit  Boundary  Procedure:  Theory  and  Appli¬ 
cations,”  AIAA  Paper  83-0123,  Reno,  NV,  January  1983. 


20 


LIST  OF  SYMBOLS 


floo 


Cm,  +  Cma 

c„ 

Cnpa 

CjV,  +  Cn^ 

^>P- 

D 


E,F,G 

H 

J 

/ 

P 


P 

Pr 

Prt 

Re 

s 

t 

u,v,w 

U,V,W 

V 


x,y,z 

^cg 


free-stream  speed  of  sound 

slope  of  the  pitching  moment  coefficient  with  angle  of  attack 
pitch-damping  moment  coefficient 
side  moment  coefficient 

slope  of  the  side  moment  coefficient  with  coning  rate 

Magnus  moment  coefficient 

pitch-damping  force  coefficient 

Magnus  force  coefficient 

projectile  diameter 

total  energy  per  unit  volume,  nondimensionalized  by 

flux  vectors  in  transformed  coordinates 

source  term  resulting  from  rotating  coordinate  frame 

Jacobian 

characteristic  length,  typically  the  projectile  diameter 
free-stream  Mach  number 

pressure,  as  used  in  thin-layer  Navier-Stokes  equations, 
nondimensionalized  by 

spin  rate  in  nonrolling  coordinate  frame,  as  used 

in  the  aerodynamic  moment  equations  and  coefficients 

Prandtl  number 

turbulent  Prandtl  number 

Reynolds  number,  OooPoo^/Moo 

viscous  flux  vector  in  transformed  coordinates 

time 

velocity  components  in  x,  y,  and  z  directions, 
nondimensionalized  by  Uqo 

Cont^av^u'iant  velocities  of  the  transformed  Navier-Stokes  equations 
free-stream  velocity  used  to  nondimensionalize  the  coning  rate, 
spin  rate,  and  the  aerodynamic  coefficients 
axial,  horizontal,  md  vertical  coordinates  with  respect  to 
the  body 

axial  location  of  projectile  center  of  gravity  with  respect  to  the 
axial  coordinate,  z 


Greek  Symbols 

a  total  angle  of  attack 

7  ratio  of  specific  heats,  as  used  in  Navier-Stokes  equations 

7  cosine  of  the  total  angle  of  attack,  as  used  in  aerodynamic 

force  and  moment  formulations 
6  sine  of  the  total  angle  of  attack 

fi  laminar  viscosity 

fit  turbulent  viscosity 
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P 

Poo 

V 

He 


transformed  coordinates  in  Navier-Stokes  equations 

density,  normalized  by  poo 

free-stream  density 

coning  rate  of  projectile 

nondimensional  coning  rate 

coning  rate  of  projectile,  nondimensionalized  by  floo  ID 
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ATTN:  Dr.  D.  W.  Hall 

994  Old  Eagle  School  Road 

Suite  1018 

Wayne,  PA  19087 

14 

Dir,  USARL 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
conunents/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-59 1 _ Date  of  Report  October  1994 _ 

2.  Date  Report  Received _ 

3.  Does  this  report,  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.) _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  doUars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 

OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City.  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  Uyre  closed,  and  mail.) 
(DO  NOT  STAPLE) 
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Postage  will  be  paid  by  addressee 

Director 

U.S.  Army  Research  Laboratory 

ATTN:  AMSRL-OP-AP-L 

Aberdeen  Proving  Ground,  MD  21005-5066 


